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SUMMARY/OVERVIEW 

The  present  research  deals  with  improving  confidence  in  the  two-scalar  mixture  fraction 
formulation  based  on  fuel  concentration  and  enthalpy  currently  applied  in  turbulent  nonpremixed 
flames.  Adding  a  third  scalar  measurement  to  the  required  fuel  Raman  and  Rayleigh  scattering  can 
reduce  uncertainty  around  the  stoichiometric  contour,  where  fuel  concentration  approaches  zero 
and  the  Rayleigh  signal  remains  relatively  constant.  Nitrogen  concentration  from  Raman  scattering 
provides  a  simultaneous,  independent,  passive  scalar,  which  has  been  used  to  modify  the  way 
functional  dependences  are  assigned  for  terms  appearing  in  the  iterative  calculation  of  mixture 
fraction.' 

TECHNICAL  DISCUSSION 

Within  the  last  year,  a  series  of  experiments  has  been  focused  on  investigating  multi-scalar 
measurements  in  turbulent  nonpremixed  flames.  A  major  goal  in  studying  such  flames  is  to 
provide  quantitative  images  of  mixture  fraction  defined  as  the  mass  fraction  of  all  atoms 
originating  from  the  fuel  stream.  This  allows  calculation  of  axial  and  radial  gradient  information, 
and  in  particular  the  scalar  dissipation  x,  which  controls  the  rate  of  molecular  mixing.  Extensive 
single-point  measurements  have  been  made,^  however  scalar  dissipation  cannot  be  extracted 
because  of  the  lack  of  spatial  gradient  data. 

One  method  for  constructing  a  conserved  scalar  suitable  for  imaging  experiments  in  reacting 
flows  has  been  through  the  simultaneous  measurement  of  temperature  (T)  and  fuel  concentration.'* 
The  conserved  scalar  j8,  is  defined  based  on  fuel  mass  fraction  (7^)  and  enthalpy  and  takes  the 
form: 

/S  =  »>  +  <:pT/Q 

where  Q  is  the  lower  heat  of  combustion  and  is  the  specific  heat  at  constant  pressure.  This  can 
be  cast  into  an  expression  for  mixture  fraction: 

_  P  -  Pair  _  ~  (^p,air\ir)lQ 

Pfuel  ”  Pair  ^FJuel  {f' p,fuel^ fuel  ~  ^p,air'^air^jQ  ^2) 


19980921  020 


This  two-scalar  approach,  which  assumes  unity  Lewis  number  and  idealized  one-step  reaction 
between  fuel  and  oxidizer,  relates  to  the  measured  signals  through: 
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where  Rm  is  the  measured  fuel  Raman  scattering  and  Ra  is  the  Rayleigh  scattering.  The  parameter 
a,  which  is  proportional  to  the  Rayleigh  cross  section,  the  mixture  molecular  weight  W,  and  the 
specific  heat  c^,  are  dependent  on  the  mixture  fraction.  Strained  counterflow  flame  calculations 
provide  appropriate  functional  forms  for  these  parameters  which  are  ineorporated  into  an  iterative 
scheme  for  determining  The  remaining  constants  Cj  and  Cj  must  be  determined  from  calibration 
experiments. 

This  approach  has  been  applied  successfully,'*'^  however  there  is  a  need  to  improve  the 
certainty  in  the  mixture  fraction  calculation  around  the  stoichiometrie  contour.  At  this  location,  the 
fuel  concentration  approaches  zero  and  the  Rayleigh  signal  remains  nearly  constant.  The  present 
work  has  examined  Raman  scattering  from  nitrogen  as  a  third  scalar  measurement  to  improve 
confidence  in  the  two-scalar  mixture  fraction  calculation. 


In  order  for  the  Nj  Raman  channel  to  provide  an  independent  conserved  scalar  with  sufficient 
signal  variation  between  regions  of  pure  air  and  pure  fuel,  experiments  were  performed  in  which 
the  fuel  stream  contained  no  nitrogen.  Assuming  no  significant  nitrogen  consumption  occurs 
during  reaction,  we  can  write  the  conserved  scalar  in  terms  of  nitrogen  mass  fraction: 

=  iW,  w 

with  the  mixture  fraction: 

where  Rrrif^  is  the  nitrogen  Raman  signal,  and  C3  is  an  additional  calibration  constant.  This 
formulation  requires  measurement  of  the  temperature  and  nitrogen  concentration,  and  thus, 
represents  an  additional  two-scalar  approach.  The  fuel  mixture  (25%  methane,  59%  argon,  and 
16%  oxygen  by  volume),  has  a  stoichiometric  mixture  fraction  ^^=0.41,  putting  the  reaction  zone 
well  inside  the  shear  layer. 

Simultaneous  planar  Rayleigh,  fuel  Raman,  and  nitrogen  Raman  images  have  been  collected  in 
experiments  using  three  cameras  and  a  single  laser.  Figure  1  shows  the  experimental  facility, 
which  is  described  in  detail  elsewhere  and  summarized  here.'  A  flashlamp-pumped  dye  laser  is 
employed  in  an  intracavity  configuration  to  generate  single-shot  energies  up  to  4.7  J  at  532  nm. 
The  beam  is  focused  into  a  sheet  over  a  6.1  mm  diameter  piloted  burner.  Scattered  light  is 
collected  on  both  sides  of  the  flame  by  low  /#  camera  lenses  oriented  perpendicular  to  the  laser 
sheet.  The  Rayleigh  scattering  and  fuel  Raman  scattering  are  colleeted  along  the  same  optical  path 
and  divided  with  a  50/50  pellicle  beam  splitter,  while  the  weaker  nitrogen  Raman  scattering  is 
collected  along  the  opposite  optical  path.  Image  intensifiers  are  lens-coupled  to  liquid-cooled  CCD 
cameras  and  isolated  with  appropriate  interference  filters. 
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In  the  past,  the  imaging  resolution  of  the  optical  setup  has  been  characterized  in  terms  of  the 
volume  associated  with  each  pixel.  For  the  present  work,  the  pixel  volume  is  49x49x500  pm^ 
where  the  largest  value  corresponds  to  the  laser  sheet  thickness,  and  the  remaining  dimensions 
describe  the  area  imaged  onto  a  single  pixel.  The  actual  spatial  resolution  is  a  more  complex 
function  of  the  optical  layout,  including  alignment,  lenses,  filters,  image  intensifiers,  and  camera 
pixel  size.  In  order  to  better  quantify  the  spatial  resolution  for  this  configuration,  simultaneous 
images  have  been  taken  of  a  uniformly  illuminated  25  pm  wire  located  at  the  focal  plane. 
Individual  camera  resolutions  are  based  upon  the  resulting  full-width  half-maximum  (FWHM) 
intensity  of  the  wire  image  from  each  camera  following  scaling,  translation,  rotation,  and  cropping. 
The  spatial  resolution  is  170  pm  on  the  Rayleigh  camera,  140  pm  on  the  nitrogen  Raman  camera, 
and  275  pm  on  the  fuel  Raman  camera.  Translating  the  wire  normal  to  the  plane  of  the  laser  sheet 
within  the  beam  thickness  (±  250  pm  from  the  focal  plane)  has  minimal  effect  (<10%)  on  the 
individual  camera  resolutions. 

Because  this  is  a  multi-camera  experiment  and  the  cameras  are  located  along  different  optical 
trains,  there  remains  the  issue  of  how  well  the  images  correlate  on  a  pixel-by-pixel  basis.  Cross¬ 
camera  spatial  resolutions  here  are  defined  based  upon  adding  combinations  of  matched  images 
from  each  camera,  and  measuring  the  resultant  full-width  half-maximum  of  the  wire.  With  optimal 
matching,  the  on-axis  spatial  resolution  is  280  pm,  about  equal  to  the  largest  single  camera 
resolution.  Factors  such  as  distortion  may  cause  degradation  away  from  the  optical  axis  to  a 
maximum  measured  400  pm.  Using  the  intensity  point  rather  than  the  FWHM  for  determining 
spatial  resolution  increases  the  reported  values  by  -75%. 

Estimates  of  the  Kolmogorov  scale  (k)  on  the  centerline  for  a  diluted  methane  flame  (3/1 
air/methane  by  volume)  give  a  value  of  k=95  pm  at  Re=20,600.'‘  Using  this  as  an  approximate 
value  for  the  turbulent  flame  in  this  work,  the  resolutions  reported  here  are  in  the  range  2-4k.  A 
study  of  an  isothermal  jet  performed  by  Namazian,  et  al.  reports  that  a  spatial  resolution  of  5k 
should  be  sufficient  for  capturing  60%  of  the  scalar  dissipation  spectrum.’  In  flames,  where  heat 
release  is  expected  to  increase  length  scales,  a  resolution  of  5k  should  be  sufficient  to  record  most 
of  the  scalar  dissipation.'* 

Figure  2  shows  line  plots  of  mixture  fraction  from  a  laminar  flame  (Re=1600)  at  a  location  15 
nozzle  diameters  downstream  (D=6. 1  mm).  Two  curves  are  shown  for  which  differ  in  the 
parameterization  of  the  mixture  fraction  dependent  terms  appearing  in  Eqn.  3  (i.e.  a,  W,  Cp).  The 
curve  marked  “No  lean  correction”  uses  flame  calculation  terms  parameterized  by  “actual”  mixture 
fraction  determined  from  the  Bilger  formula.^ 

Figure  3  shows  that  this  approach  incorporates  a  departure  of  from  the  more  rigorous 
formulation,  as  indicated  by  strained  laminar  flame  calculations  (strain  rate=100  s  ’).  Deviation 
from  one-step  chemistry  (i.e.  loss  of  parent  fuel  to  intermediate  species)  is  compensated  in  the 
curve  marked  by  using  a  weighting  term  involving  reactivity.®  The  curve  for  exhibits  little 
deviation  from  the  actual  mixture  fraction.  By  assigning  functional  dependences  based  on  the 
predicted  ^  from  flame  calculations,  the  curve  shown  in  Fig.  2  is  obtained  showing 
improved  agreement  with  ^*^2.  The  flame  computations  are  insensitive  to  variations  in  the  strain  rate 


over  the  range  10-200  s  ',  which  is  expected  to  be  representative  of  the  scalar  dissipation  values 
measured  in  the  turbulent  flame,  based  on  existing  data  in  similar  flames.*  '" 

Applying  this  technique  to  single-pulse  imaging  in  turbulent  nonpremixed  flames  provides 
similar  results,  although  the  N2  Raman  signal  is  affected  by  noise,  especially  in  regions  of  high 
mixture  fraction  where  there  is  little  nitrogen.  Figure  4  shows  images  taken  25D  downstream  from 
a  Re=  15,000  turbulent  flame;  the  Raman  images  have  been  contour  smoothed."  Qualitatively,  the 
scalar  dissipation  fields,  %  (defined  as  =  2^  - V^,  where  2)  is  the  diffusivity),  appear 

similar,  revealing  the  same  main  structural  features,  with  significant  scalar  dissipation  apparent 
along  the  edge  of  the  main  jet.  The  position  of  the  stoichiometric  mixture  fraction  contour  is 
highlighted  (black  lines)  in  these  images. 

In  summary,  recent  experiments  measuring  nitrogen-temperature  mixture  fraction  have  served 
as  a  guide  for  correcting  the  fuel-temperature  mixture  fraction  for  values  around  and  lean  of 
stoichiometric.  It  has  been  shown  that  parameterizing  specific  heat,  molecular  weight,  and 
Rayleigh  cross  section  as  a  function  of  predicted  from  counterflow  flame  calculations,  rather 
than  the  actual  mixture  fraction,  improves  the  performance  of  this  two-scalar  approach.  Under 
turbulent  conditions  (i?e=  15,000),  the  two  approaches  reveal  differences  close  to  the  centerline, 
most  likely  a  result  of  noise  limitations  of  the  nitrogen  Raman  signal.  This  work  increases 
confidence  in  employing  ^  for  mixture  fraction  determination,  which  remains  the  most  attractive 
approach  because  of  its  superior  signal-to-noise  characteristics.  Further  details  of  the  results  are 
documented  in  the  paper  accepted  to  the  Twenty-Seventh  Symposium  on  Combustion,  which  will 
be  published  later  this  year.' 
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CH4  Raman  qqP  cameras 

Figure  1.  Schematic  of  the  three-scalar  intracavity 
turbulent  flame  imaging  experiment. 


Radial  Location  (r/D) 

Figure  2.  Radial  variation  of  mixture  fraction  calculated 
from  fuel-temperature  (solid  line)  and  nitrogen- 
temperature  (short  dashes)  two-scalar  approaches  15D 
downstream  in  a  laminar  flame  (Re=1600). 
overpredicts  the  fuel-temperature  mixture  fraction  in 
regions  around  and  lean  of  stoichiometric  (long  dashes). 


Actual  Mixture  Fraction 

Figure  3.  Mixture  fraction  calculated  from  strained 
laminar  flame  calculations  (100  s‘^)  using  the  fiiel- 
temperature  (solid  line)  and  nitrogen-temperature  (short 
dashes)  two-scalar  approaches  plotted  against  mixture 
fraction  calculated  using  the  formula  proposed  by  Bilger 
[7].  The  effect  of  fuel  correction  on  is  shown  (long 
dashes). 
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Figure  4.  Instantaneous  (a)  Rayleigh,  (b)  CH4  Raman, 
(c)  N2  Raman  images  of  the  turbulent  flame  taken  25 
jet  diameters  (D=6.1  mm)  downstream.  Also  shown  are 
the  computed  mixture  fraction  images  (d)  (e)  ^^2, 

(f)  temperature,  and  scalar  dissipation  images  (g) 
and  (h)  Each  image  size  is  4.3D  x  ID. 
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Abstract 


Raman  scattering  from  nitrogen  in  laminar  (Rc=1600)  and  turbulent  (/?e=  15,000)  nonpremixed 
flames  of  diluted  methane  has  been  investigated  as  a  means  of  complementing  the  two-scalar 

mixture  fraction  computation  based  on  measurements  of  fuel  mass  fraction  and  temperature 

Using  a  diluent  consisting  of  argon  and  oxygen  with  an  overall  volumetric  dilution  ratio  of  3/1 
(diluent/fuel)  allows  sufficient  variation  in  the  measured  nitrogen  concentration  for  its  use  as  a 
passive  scalar.  The  present  experimental  setup  requires  only  a  single  laser  (532  nm)  in  a  high- 
power  intracavity  configuration.  Mixture  fraction  profiles  calculated  using  independent  fuel- 
temperature  and  nitrogen-temperature  two-scalar  approaches  show  excellent  agreement  in  the 
laminar  flame.  For  the  turbulent  flame,  a  100  single-shot  average  of  mixture  fraction  shows 
reasonable  agreement  between  the  two  approaches.  Discrepancies  are  most  evident  in  single-shot 

images  in  regions  of  large  mixture  fraction  (^>^3)  where  the  nitrogen  Raman  signal  is  noise- 

dominated.  The  location  of  the  stoichiometric  contour  is  consistently  determined  by  both  two- 
scalar  approaches.  A  simple  modification  of  the  functional  dependence  for  parameters  appearing  in 
the  fuel-temperature  mixture  fraction  formulation  is  shown  to  correct  deficiencies  in  the  approach 

predicted  by  laminar  flame  calculations,  most  notably  for  values  of  around  stoichiometric  and 
into  lean  regimes. 


Introduction 


Planar  imaging  of  turbulent  flames  using  laser  diagnostic  techniques  such  as  Lorenz-Mie  [1], 
Rayleigh  [2],  laser-induced  fluorescence  [3],  and  Raman  scattering  [4],  have  long  been  directed 
towards  improved  characterization  of  the  flow  structures  and  flame  front  location.  One  major  goal 

in  studying  turbulent  flames  is  to  provide  quantitative  images  of  mixture  fraction  defined  as  the 
mass  fraction  of  all  atoms  originating  from  the  fuel  stream.  This  allows  calculation  of  axial  and 
radial  gradient  information,  and  in  particular  the  scalar  dissipation  which  controls  the  rate  of 

molecular  mixing.  Experimental  data  are  essential  for  modelers  to  construct  joint  probability 
density  functions  conditional  on  these  parameters.  Extensive  single-point  measurements  have  been 
made  [5],  however  scalar  dissipation  cannot  be  extracted  because  of  the  lack  of  spatial  gradient 
data.  Line-Raman  measurements  [6,7]  are  less  restrictive,  allowing  the  computation  of  one¬ 
dimensional  scalar  dissipation. 

One  method  for  constructing  a  conserved  scalar  suitable  for  reacting  flows  has  been  through 
the  simultaneous  measurement  of  temperature  (T)  and  fuel  concentration  [8].  The  conserved  scalar 

P,  is  defined  based  on  fuel  mass  fraction  (Y^)  and  enthalpy  and  takes  the  form: 

P  =  Yp  +  CpT/Q 

where  Q  is  the  lower  heat  of  combustion  and  is  the  specific  heat  at  constant  pressure.  This  can 
be  cast  into  an  expression  for  mixture  fraction: 

eFT  _  P  ~  Pair  _  ~  '^air)/Q 

P fuel  ~  Pair  ^FJuel  fuel  ~  ^air  ) / 2 


This  two-scalar  approach,  which  assumes  unity  Lewis  number  and  idealized  one-step  reaction 
between  fuel  and  oxidizer,  relates  to  the  measured  signals  through: 
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where  Rm  is  the  measured  fuel  Raman  scattering  and  Ra  is  the  Rayleigh  scattering.  The  parameter 
a,  which  is  proportional  to  the  Rayleigh  cross  section,  the  mixture  molecular  weight  W,  and  the 
specific  heat  c^,  are  dependent  on  the  mixture  fraction.  Strained  counterflow  flame  calculations 
provide  appropriate  functional  forms  for  these  parameters  which  are  incorporated  into  an  iterative 


scheme  for  determining  The  remaining  constants  Cj  and  Cj  must  be  determined  from  calibration 
experiments. 

This  approach  has  been  applied  successfully  [3,  9-14],  however  there  is  a  need  to  improve  the 
certainty  in  the  mixture  fraction  calculation  around  the  stoichiometric  contour.  At  this  location,  the 
fuel  concentration  approaches  zero  and  the  Rayleigh  signal  remains  nearly  constant.  The 
acquisition  of  an  appropriate  third  scalar  with  continuously  varying  signal  could  improve 
confidence  in  the  two-scalar  mixture  fraction  calculation  in  this  regime  by  providing  independent 
quantitative  confirmation  of  the  measurements.  The  present  work  examines  Raman  scattering  from 
nitrogen  as  a  third  scalar  measurement  in  turbulent  nonpremixed  flames. 

In  this  paper,  simultaneous  planar  Rayleigh,  fuel  Raman,  and  nitrogen  Raman  images  have 
been  collected  in  experiments  using  three  cameras  and  a  single  laser.  Current  detector  technology 
combined  with  high-energy  (~5  J/pulse)  laser  configurations  makes  Ramanography  feasible  for 
molecules  sucb  as  nitrogen,  despite  its  small  scattering  cross-section  [15].  In  order  for  the  N2 
Raman  channel  to  provide  an  independent  passive  conserved  scalar  with  sufficient  signal  variation 
between  regions  of  pure  air  and  pure  fuel,  experiments  were  performed  in  which  the  fuel  stream 
contained  no  nitrogen.  Assuming  no  significant  nitrogen  consumption  occurs  during  reaction,  we 
can  write  the  conserved  scalar  in  terms  of  nitrogen  mass  fraction: 


with  the  mixture  fraction: 
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where  Rntj^^  is  the  nitrogen  Raman  signal,  and  C3  is  an  additional  calibration  constant.  This 

formulation  requires  measurement  of  the  temperature  and  nitrogen  concentration,  and  thus, 
represents  an  additional  two-scalar  approach. 

Experiment 

The  experimental  facility  used  in  the  present  work  to  measure  Rayleigh,  fuel  Raman,  and 
nitrogen  Raman  scattering  is  shown  schematically  in  Fig.  1.  A  single  flashlamp-pumped  dye  laser 
(Candela  LFDL-20;  Pyrromethene  546  dye,  2.0xl0'^M  in  methanol)  is  employed  in  an  intracavity 
configuration  to  generate  single-shot  energies  up  to  4.7  J  at  532  nm.  The  beam  is  focused  into  a 
sheet  by  a  30-cm  focal-length  cylindrical  lens,  and  an  8  mm  wide  slit  inserted  into  the  cavity 
provides  a  beam  waist  of  500  pm.  Beam  waist  measurements  are  performed  by  imaging  the  laser 
sheet  from  above  using  Lorenz-Mie  scattering  from  NaCl  particles  seeded  into  air  with  a  nebulizer 
(TSI  Model  9306). 

Scattered  light  is  collected  on  both  sides  of  the  flame  by  low  /#  camera  lenses  oriented 

perpendicular  to  the  laser  sheet.  The  Rayleigh  scattering  and  fuel  Raman  scattering  are  collected 
along  the  same  optical  path  and  divided  with  a  50/50  pellicle  beam  splitter,  while  the  weaker 
nitrogen  Raman  scattering  is  collected  along  the  opposite  optical  path.  Image  intensifiers  are  lens- 
coupled  to  liquid-cooled  CCD  cameras  and  isolated  with  appropriate  10  nm  bandwidth  interference 
filters.  The  filters  and  cameras  used  are:  532  nm  (center  wavelength)  for  the  Rayleigh  scattering 
(Photometries  Starl  CCD),  630  nm  for  the  Stokes-shifted  Raman  scattering  from  methane 
(Photometries  CH250),  and  610  nm  for  nitrogen  Raman  scattering  (Princeton  Instruments 
TE/CCD-512).  The  image  intensifiers  are  gated  for  the  same  2  |is  period  to  bracket  the  laser  pulse 
and  minimize  background  interferences  from  flame  luminosity.  The  experiment  is  controlled  from 
a  single  computer  which  handles  the  subsequent  image  transfer  and  storage.  Raw  images  are 
corrected  for  background  and  throughput  as  described  in  Ref.  16. 
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In  the  past,  the  imaging  resolution  of  the  optical  setup  has  been  characterized  in  terms  of  the 
volume  associated  with  each  pixel.  For  the  present  work,  the  pixel  volume  is  49x49x500  pm\ 
where  the  largest  value  corresponds  to  the  laser  sheet  thickness,  and  the  remaining  dimensions 
describe  the  area  imaged  onto  a  single  pixel.  The  actual  spatial  resolution  is  a  more  complex 
function  of  the  optical  layout,  including  alignment,  lenses,  filters,  image  intensifiers,  and  camera 
pixel  size.  In  order  to  better  quantify  the  spatial  resolution  for  this  configuration,  simultaneous 
images  have  been  taken  of  a  uniformly  illuminated  25  pm  wire  located  at  the  focal  plane. 
Individual  camera  resolutions  are  based  upon  the  resulting  full-width  half-maximum  (FWHM) 
intensity  of  the  wire  image  from  each  camera  following  scaling,  translation,  rotation,  and  cropping. 
The  spatial  resolution  is  170  pm  on  the  Rayleigh  camera,  140  pm  on  the  nitrogen  Raman  camera, 
and  275  pm  on  the  fuel  Raman  camera.  Translating  the  wire  normal  to  the  plane  of  the  laser  sheet 
within  the  beam  thickness  (±  250  pm  from  the  focal  plane)  has  minimal  effect  (<10%)  on  the 
individual  camera  resolutions. 

Because  this  is  a  multi-camera  experiment  and  the  cameras  are  located  along  different  optical 
trains,  there  remains  the  issue  of  how  well  the  images  correlate  on  a  pixel-by-pixel  basis.  Cross¬ 
camera  spatial  resolutions  here  are  defined  based  upon  adding  combinations  of  matched  images 
from  each  camera,  and  measuring  the  resultant  full-width  half-maximum  of  the  wire.  With  optimal 
matching,  the  on-axis  spatial  resolution  is  280  pm,  about  equal  to  the  largest  single  camera 
resolution.  Factors  such  as  distortion  may  cause  degradation  away  from  the  optical  axis  to  a 
maximum  measured  400  pm.  Using  the  Me^  intensity  point  rather  than  the  FWHM  for  determining 
spatial  resolution  increases  the  reported  values  by  -75%. 

Estimates  of  the  Kolmogorov  scale  (k)  on  the  centerline  for  a  diluted  methane  flame  (3/1 

air/methane  by  volume)  give  a  value  of  k=95  pm  at  Re=20,600  [9].  Using  this  as  an  approximate 
value  for  the  turbulent  flame  in  this  work,  the  resolutions  reported  here  are  in  the  range  2-4k.  A 
study  of  an  isothermal  jet  performed  by  Namazian,  et  al.  [17]  reports  that  a  spatial  resolution  of  5k 
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should  be  sufficient  for  capturing  60%  of  the  scalar  dissipation  spectrum.  In  flames,  where  heat 
release  is  expected  to  increase  length  scales,  a  resolution  of  5k  should  be  sufficient  to  record  most 

of  the  scalar  dissipation  [9]. 

The  burner  consists  of  a  6.1  mm  diameter  nozzle  surrounded  by  a  14  mm  pilot  flame  region, 
which  prevents  blowoff  of  the  turbulent  flames  and  allows  examination  of  a  large  range  of 
Reynolds  numbers.  The  premixed  pilot  flame  is  fueled  by  stoichiometric  acetylene/hydrogen/air  in 
proportions  such  that  the  mixture  carbon/hydrogen  ratio  matches  that  of  the  main  jet,  which  allows 
the  flame  to  be  modeled  as  a  two-stream  mixing  problem  [18].  A  filtered  ambient  air  coflow 
(~7  m/s)  prevents  entrainment  of  dust  into  the  flame  that  would  interfere  with  the  Rayleigh 
measurement.  The  fuel  used  in  this  work  is  methane  mixed  with  argon  and  oxygen  to  provide  a 
3/1  dilution  ratio  by  volume  and  a  diluent  oxygen  content  to  match  that  of  air.  For  this  fuel 

composition,  the  stoichiometric  mixture  fraction  is  =  0.41  which  puts  the  reaction  zone  well 
inside  the  shear  layer. 

Dilution  suppresses  much  of  the  soot  and  formation  of  soot  precursors  that  can  affect  the 
nitrogen  Raman  signal  [19].  In  the  laminar  flame,  Cj  fluorescence  contributes  about  15%  to  the 
Raman  signal  on  the  rich  side  of  the  flame  front  when  isolated  with  the  10  nm  bandwidth 
interference  filter.  Using  a  narrower  3  nm  bandwidth  filter  centered  at  a  wavelength  of  607.3  nm 
virtually  eliminates  this  interference  at  a  loss  of  half  of  the  Raman  signal.  The  3  nm  filter  is  used 
for  the  laminar  flame  results  presented  in  the  next  section.  The  broader  10  nm  bandwidth  filter  is 
used  for  the  turbulent  flame  imaging  to  provide  maximum  SNR,  as  the  interference  is  expected  to 
be  less  significant  than  under  laminar  flame  conditions. 

The  N2  Raman  images  are  corrected  for  temperature  dependence  of  the  scattering  cross  section 
and  the  overlap  of  the  Raman  spectral  profile  with  the  bandpass  of  the  interference  filter.  Simple 
convolution  of  the  temperature-dependent  Raman  spectmm  (computed  using  the  Raman  code 
developed  by  Hassel  [20])  with  the  zero-degree-incidence  filter  profile  gives  a  correction  factor  that 
does  not  adequately  predict  the  signal  variation  with  temperature.  This  effect  is  a  consequence  of 
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the  placement  of  the  interference  filter  relative  to  the  collection  lens.  A  geometric  analysis  reveals 
that  off-axis  rays  entering  the  filter  create  a  significant  shift  in  the  effective  spectral  profile  towards 
shorter  wavelengths.  A  series  of  calibration  experiments,  in  conjunction  with  the  geometric 
analysis,  is  used  to  provide  the  temperature  correction.  The  magnitude  of  the  Nj  Raman  correction 

is  as  high  as  50%  at  the  peak  flame  temperatures  (13^-2200  K)  for  the  10  nm  bandwidth 

interference  filter.  This  temperature  dependence  is  partly  beneficial  in  that  ample  N2  Raman  signal 
is  available  even  at  higher  temperatures.  A  similar  analysis  applied  to  the  methane  Raman  signal 
indicates  no  off-axis  ray  correction  is  necessary. 

Calibration  images  are  taken  before  and  after  each  experimental  sequence  in  flows  of  known 
uniform  concentration  (helium,  methane,  and  air).  Additionally,  a  series  of  measurements 
performed  in  a  non-reacting  jet  of  fuel  issuing  into  ambient  air  (Tj=T2=298  K)  confirms  the 
precision  of  the  temperature  images,  with  maximum  fluctuations  about  the  mean  of  ±  10  K  and  a 
signal-noise  ratio  (SNR)  of  70. 

Results  and  Discussion 

Images  have  been  obtained  from  a  laminar  flame  (Re=1600)  taken  15  nozzle  diameters  (D) 
downstream  of  the  jet  exit,  and  from  a  turbulent  flame  (Re=  15,000)  at  a  distance  25D  downstream. 
Data  from  the  laminar  flame  represent  a  50-shot  average,  which  gives  a  Rayleigh  SNR  of  250  in 
the  ambient  air  region.  Fuel  Raman  SNR  is  30  in  the  fuel  core,  and  the  nitrogen  Raman  has  SNR 
values  of  20  in  the  ambient  air.  A  contour  smoothing  technique  [21]  has  been  applied  to  both 
Raman  channels,  which  improves  the  fuel  Raman  SNR  by  a  factor  of  2  and  the  nitrogen  Raman  by 
more  than  a  factor  of  10. 

Figure  2  shows  the  scattered  signals  along  a  line  (one  pixel  height)  as  a  function  of  the 
normalized  radius.  The  flattening  of  the  N2  Raman  profile  around  r/D  =  0.75  is  attributable  to  tbe 
temperature  sensitivity  of  the  signal;  increased  spectral  overlap  with  the  interference  filter  occurs  at 
higher  temperatures  for  this  experimental  setup.  Figure  3  shows  mixture  fraction  as  a  function  of 
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radius  in  the  laminar  flame  computed  using  fuel-temperature  and  nitrogen-temperature 
two-scalar  approaches.  There  is  good  agreement  between  the  two  profiles,  and  both  predict  a 
centerline  mixture  fraction  just  over  ^-0.9,  and  the  radial  location  of  the  stoichiometric  contour 
coincides. 

Counterflow  flame  calculations  (strain  rate=100  s  ')  of  this  system  demonstrate  the  departure  of 
the  two-scalar  formulation  from  the  “actual”  mixture  fraction  (Fig.  4)  calculated  using  the  Bilger 
formula  [22]: 

IZcjWc  +  +(Zo.»  -  ZolM 

^^C,fuel  I^C  T ^Hyfuel /^H  {^0,air  ~  ^0,fuel 

where  Zj  and  Wj  are  the  atomic  mass  fractions  and  weights,  respectively,  for  carbon  (C),  hydrogen 
(H),  and  oxygen  (O).  Deviation  from  one-step  chemistry  (i.e.  loss  of  parent  fuel  to  intermediate 
species)  is  compensated  in  calculating  fuel  concentration  by  the  expression  Rnif  =  (1  +  Cp 

O^),  where  Rm-p,^  is  the  uncompensated  value,  and  the  reactivity  is  defined  as  O  =  (T  -  T3„bientV(Tad 
-Ta„bie„.)[12]-  Cp  is  a  weighting  coefficient  (Cp  =  0.8)  and,  as  illustrated  in  the  figure  by  the  curve 
marked  corrects  on  the  rich  side  of  stoichiometric.  The  slight  plateau  in  mixture  fraction 

shown  in  Figure  3  coincides  with  the  deviation  occurring  near  ^"^=0.5  in  Fig.  4,  and  spatially 
coincides  with  the  peak  temperature.  It  is  also  predicted  from  the  flame  calculations  that  the 
mixture  fraction  based  upon  nitrogen  mass  fraction  should  exhibit  little  deviation  from  To 

correct  on  the  lean  side,  Cp(^),  o(^),  and  W(^)  are  described  as  functions  of  the  predicted  two- 

scalar  mixture  fraction  rather  than  the  actual  mixture  fraction  from  the  flame  calculations.  The 

latter  approach  gives  the  third  curve  in  Fig.  3  (denoted  F-T  no  lean  correction),  which  overpredicts 
the  mixture  fraction  —  in  agreement  with  the  trends  indicated  in  Fig.  4.  Without  this  simple 
modification,  the  location  of  the  stoichiometric  contour  is  in  error  by  O.ID  towards  the  lean  side. 
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The  temperature  correction  of  the  N2  Raman  signal  causes  a  minor  dependence  of  on  the  value 

of  because  the  latter  is  used  to  adjust  the  Rayleigh  cross  section  parameter,  a(^).  This  effect  is 

small,  however,  and  the  magnitude  of  the  temperature  correction  is  reduced  by  20%  in  the 
turbulent  flame  measurement  because  of  the  larger  bandwidth  interference  filter.  The  flame 
computations  are  insensitive  to  variations  in  the  strain  rate  over  the  range  10-200  s’’,  which  is 
expected  to  be  representative  of  the  scalar  dissipation  values  measured  in  the  turbulent  flame,  based 
on  existing  data  in  similar  flames  [12, 14, 21]. 

For  the  turbulent  flame,  100  instantaneous  images  were  processed  and  then  averaged.  Figure  5 
shows  the  Rayleigh  and  Raman  profiles.  Comparison  of  the  computed  mixture  fractions  in  Fig.  6 
shows  good  agreement  between  the  two  approaches,  and  the  location  of  the  stoichiometric  contour 

again  coincides.  For  values  ^>^5,  there  is  some  discrepancy  between  the  profiles  of  and 

and  peak  centerline  mixture  fractions  are  ^|^o=0.75  versus  ^|l'o=0.80.  These  discrepancies  may 

reflect  the  limitations  in  ^  caused  by  the  deviation  from  the  actual  mixture  fraction  indicated  in 

Fig.  4;  however,  the  maximum  discrepancy  between  the  two  profiles  has  a  width  of  only  0.2£). 
Again,  failure  to  parameterize  the  variables  in  the  two-scalar  mixture  fraction  with  respect  to 

predicted  would  have  caused  a  larger  discrepancy  in  the  mixture  fraction  obtained  with  the 
different  two-scalar  techniques. 

In  order  to  examine  the  performance  of  the  different  mixture  fraction  formulations  on  a  single¬ 
shot  basis  in  the  turbulent  flame,  one  set  of  collected  signals  is  shown  in  Fig.  7  and  the  derived 

mixture  fractions  are  shown  in  Fig.  8.  Line  plots  of  mixture  fraction  reveal  that  and  profiles 

generally  agree  around  stoichiometric,  although  noise  limitations  are  evident  in  the  nitrogen  Raman 
signal.  In  the  turbulent  flame,  the  Raman  signal  has  a  SNR  of  6  before  smoothing;  this 
increases  to  30  after  contour  smoothing.  Instantaneous  images  of  the  turbulent  flame  are  shown  in 
Fig.  9,  with  an  image  size  of  4.3D  x  ID.  Radial  data  for  Figs.  7-8  are  taken  along  the  center  of 


8 


each  image.  Additionally,  scalar  dissipation  (two-dimensions)  is  computed  based  upon  each 
mixture  fraction  formulation,  where  the  scalar  dissipation  %,  is  defined  as  ;i;  =  2.0  V<^. 

Temperature  variations  are  accounted  for  in  the  diffusivity  by  the  relationship  (T/Tq)' 

The  image  should  be  regarded  as  semi-quantitative  due  to  the  modest  SNR  of  the  nitrogen 

Raman  signal.  Qualitatively,  the  images  appear  similar,  revealing  the  same  main  structural 
features,  with  peak  scalar  dissipation  apparent  along  the  edge  of  the  main  jet.  The  position  of  the 
stoichiometric  mixture  fraction  contour  is  highlighted  (black  lines)  in  these  images. 


Conclusions 

Two  simultaneous  methods  of  imaging  mixture  fraction  have  been  investigated  using  planar 
Rayleigh,  fuel  Raman,  and  nitrogen  Raman  scattering  in  nonpremixed  flames  using  argon/oxygen 
diluted  methane  as  a  fuel.  The  different  approaches  are  in  good  agreement  when  applied  to  a 
laminar  flame.  Here,  the  nitrogen-temperature  mixture  fraction  served  as  a  guide  for  correcting  the 
fuel-temperature  mixture  fraction  for  values  around  and  lean  stoichiometric.  It  has  been  shown  that 

parameterizing  specific  heat,  molecular  weight,  and  Rayleigh  cross  section  as  a  function  of 

predicted  from  counterflow  flame  calculations,  rather  than  the  actual  mixture  fraction,  dramatically 
improves  the  performance  of  this  two-scalar  approach.  Under  turbulent  conditions  (/?e=  15,000), 
the  two  approaches  reveal  differences  close  to  the  centerline,  most  likely  a  result  of  noise 
limitations  of  the  nitrogen  Raman  signal.  One  of  the  shortcomings  of  the  two-scalar  approach 
using  fuel  concentration  and  Rayleigh  scattering  is  the  lack  of  a  strongly  varying  scalar  quantity 
around  the  stoichiometric  contour.  This  work  has  shown  that  Stokes-shifted  nitrogen  Raman 
scattering  is  a  viable  candidate  for  a  third  scalar  quantity  to  supplement  the  usual  two-scalar 
approach  based  on  measurements  of  fuel  mass  fraction  and  temperature.  By  eliminating  nitrogen 
from  the  fuel  stream  diluent,  the  nitrogen  mass  fraction  provides  sufficient  variation  to  serve  as  a 
passive  scalar. 
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Although  modest  signal-noise  of  the  nitrogen  Raman  make  more  detailed  quantitative  and 
statistical  comparisons  premature  at  this  point,  the  current  evidence  indicates  that  this  technique 
merits  continued  investigation.  An  added  benefit  of  this  experimental  setup  is  that  only  a  single 
laser  is  required.  The  ultimate  goal  of  these  efforts  is  to  combine  three  measurements  into  a 
comprehensive  mixture  fraction  calculation  which  offers  significant  improvement  over  two-scalar 
techniques,  while  maintaining  the  ability  to  compute  gradient  information  required  to  accurately 
determine  scalar  dissipation. 
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Figure  Captions 


1 .  Schematic  of  the  three-scalar  turbulent  flame  imaging  experiment. 

2.  Line  plot  (one  pixel  height)  showing  radial  variation  of  Rayleigh,  fuel  Raman,  and  nitrogen 
Raman  signals  in  the  laminar  flame  (Rc=1600)  15  jet  diameters  (D=6A  mm)  downstream  of  the  jet 
exit.  The  data  represent  a  50-shot  average. 

3 .  Radial  variation  of  mixture  fraction  calculated  from  fuel-temperature  (solid  line)  and  nitrogen- 
temperature  (short  dashes)  two-scalar  approaches  15D  downstream  in  the  laminar  flame.  Previous 

functional  dependences  of  parameters  in  overpredict  the  fuel-temperature  mixture  fraction  in 
regions  around  and  lean  of  stoichiometric  (long  dashes). 

4 .  Mixture  fraction  calculated  from  strained  laminar  flame  calculations  (100s  ')  using  the  fuel- 
temperature  (solid  line)  and  nitrogen-temperature  (short  dashes)  two-scalar  approaches  plotted 
against  mixture  fraction  calculated  using  the  formula  proposed  by  Bilger  [22].  The  effect  of  fuel 

correction  on  is  also  shown  (long  dashes). 

5 .  Measured  radial  Raman  and  Rayleigh  intensities  for  a  100-shot  average  of  the  turbulent  flame 
(Rc=  15,000)  along  a  line  25D  downstream. 

6 .  Radial  variation  of  mixture  fraction  calculated  from  fuel-temperature  (solid  line)  and  nitrogen- 
temperature  (short  dashes)  two-scalar  approaches  25D  downstream  in  the  turbulent  flame  (100- 
shot  average). 

7 .  Measured  radial  Raman  and  Rayleigh  intensities  for  a  single  shot  of  the  turbulent  flame  along  a 
line  25D  downstream. 

8 .  Radial  variation  of  mixture  fraction  calculated  from  fuel-temperature  (solid  line)  and  nitrogen- 
temperature  (short  dashes)  two-scalar  approaches  25D  downstream  in  the  turbulent  flame  (single 
shot). 

9 .  Instantaneous  (a)  Rayleigh,  (b)  CH4  Raman,  (c)  N2  Raman  images  of  the  turbulent  flame. 

Also  shown  are  the  computed  mixture  fraction  images  (d)  (e)  (f)  temperature,  and  scalar 

dissipation  images  (g)  and  (h)  Image  size  is  4.3D  x  ID;  data  for  line  plots  of  Figs.  7-8 
are  t^en  along  the  center  of  each  image. 
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Simultaneous  two-dimensional  Rayleigh  and  fuel  Raman  images  have  been  collected  m  air- 
diluted  methane  and  hydrogen  jet  diffusion  flames.  Temperature,  fuel  mass  fraction  and  mixture 
fraction  images  are  derived  by  a  two-scalar  approach  based  on  one-step  chemistry  and  equal 
species  diffusivities.  This  enables  calculation  of  two  components  of  the  scalar  dissipation  rate  x- 
The  inherently  weak  Raman  signal  has  been  maximised  by  intra-cavity  measurernents,  using  a 
flashlamp-pumped  dye  laser.  In  addition,  the  Raman  signal-to-noise  ratio  is  drastically 
improved  by  a  novel  contour-aligned  smoothing  technique  which  exploits  the  high  correlation 
between  the  Rayleigh  and  Raman  signals. 

Quantitative  measurements  of  scalar  dissipation  are  presented,  including  probability  density 
functions  for  components  of  x-  Profiles  of  mean  and  rms  mixture  fraction  show  the  usual 
features  already  documented  in  other  published  results  for  this  type  of  flame.  Probability 
density  functions  of  ^  are  close  to  Gaussian  on  the  axis,  and  tend  to  bimodal  at  the  edge  of  the 
flame.  Results  for  the  CH4  flames  indicate  that  the  mean  of  x  shows  little  change  with  Reynolds 
number.  In  the  H2  flame,  mean  values  for  the  axial  and  radial  components  of  the  scalar 
dissipation  rate,  x»  ^ire  nearly  the  same,  indicating  a  more  isotropic  structure  than  in  the  CH4 
flames.  For  both  fuels,  the  pdf  of  ln(x)  on  the  axis  is  more  peaky  than  a  lognormal  distribution 
and  somewhat  skewed.  The  profiles  of  (xW)  show  a  nonlinear  dependence  on  mixture  l^raction 
and  have  no  clear  resemblance  to  the  skewed,  monomodal  shapes  seen  in  cold  flows.  In  the  Hj 
flame  there  is  a  strong  correlation  between  instantaneous,  local  values  of  scalar  dissipation  and 
the  departure  from  equilibrium,  as  measured  by  temperature  depression. 

Keywords:  Laser  diagnostics;  turbulent  nonpremixed  flames;  scalar  dissipation 


’’’Corresponding  author.  E-mail:  bilger(§tiny.me.su. oz.au 


141 


142 


S.  H.  STARNER  et  at. 


1.  INTRODUCTION 

In  an  overview  of  recent  development  of  techniques  (Frank  et  al.,  1995, 
Kelmen  et  al.,  1995;  Starner  et  al.,  1994b)  for  measuring  mixture  fraction 
and  its  dissipation  in  turbulent  flames,  (Bilger,  1993a)  argues  that  such 
experiments  are  needed  not  only  to  test  the  validity  of  flamelet  models,  but 
also  to  aid  in  the  development  of  the  Conditional  Moment  Closure 
approach  ( Klimenko,  1990;  Bilger,  1993b),  where  information  on  the  mean 
scalar  dissipation  (xlO.  conditional  on  the  mixture  fraction,  is  vital. 
Bilger  (1993a)  also  stresses  that  only  by  the  use  of  a  conserved  scalar  such  as 
^  can  turbulent  flame  structure  and  statistics  be  directly  related  to  non- 
reacting  turbulent  flows. 

The  two-scalar  theory  (which  uses  the  simplifying  assumptions  of  one- 
step  chemistry  and  equal  species  diffusivities)  on  which  the  present  work  is 
based  (Bilger,  1993c)  has  been  validated  by  single-point  Raman/Rayleigh 
measurements  (Starner  et  al.,  1992).  There  are  several  possible  combinations 
of  scalar  pairs:  in  early  experiments,  a  fluorescing  fuel  (acetaldehyde)  was 
used,  which  in  combination  with  Rayleigh  scattering  gave  a  pair  of  strong 
signals  with  a  high  signal-to-noise  ratio  (SNR)  ( Starner  et  al.,  1994b,  Frank 
et  al.,  1995).  This  enabled  good  scalar  dissipation  data  to  be  derived,  but 
loss  of  parent  fuel  to  intermediates  was  found  to  cause  too  large  a  departure 
from  the  one-step  chemistry  assumption  implicit  in  the  two-scalar  method. 
Not  only  does  fuel  pyrolysis  contribute  to  this  loss,  but  also  the  inherently 
high  CO  production  in  oxygen-containing  fuels  such  as  acetaldehyde, 
acetone  and  alcohols  make  these  fuels  less  compatible  with  the  two-scalar 
method. 

An  alternative  is  to  record  directly  the  spontaneous  Raman  signal  from  an 
alkane  fuel  such  as  CH4.  This  has  been  implemented  successfully  by  a 
modified  line  imaging  technique  (Starner  et  al.,  1994a).  However,  for  2D 
imaging  the  weak  Raman  signal  has  until  recently  proved  too  noisy  for  the 
purpose  of  obtaining  x  quantitatively  (Kelman  et  al.,  1995).  To  overcome 
this  problem,  an  intracavity  layout  is  now  in  use  to  increase  the  laser  sheet 
energy,  but  even  so,  the  SNR  in  the  raw  Raman  image  is  rarely  above  eight. 
This  is  far  from  adequate  for  producing  even  qualitative  scalar  dissipation 
profiles;  the  Raman  SNR  needs  to  be  above  30  for  this  purpose.  Some 
improvement  can  be  obtained  with  conventional  smoothing  (indiscrimi¬ 
nately  averaging  over  a  fixed  size  sub-domain,  ‘  top-hat  or  Gaussian)  but 
the  smoothing  required  to  produce  realistic-looking  scalar  dissipation 
profiles  is  so  extensive  that  a  large  part  of  the  gradient  information  is  lost; 
the  dissipation  would  be  at  best  qualitative.  However,  we  have  exploited  the 
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high  correlation  between  the  scalars  used  here  (Rayleigh  and  fuel  Raman)  to 
align  the  smoothing  sub-domain  in  the  Raman  image  along  the  contours  of 
the  Rayleigh  image.  It  is  shown  (Starner  et  al.,  1995)  that  up  to  a  factor  of 
ten  improvement  in  Raman  SNR  is  obtainable  by  this  method.  With  the 
present  combination  of  intracavity  layout  and  contour-aligned  smoothing, 
we  are  able  here  to  obtain  adequate  SNR  even  with  the  very  weak  Raman 
signal  of  the  H2  flame.  The  H2  flame  has  the  advantages  that  there  is  no  soot 
problem  and  virtually  no  fuel  loss  to  intermediates;  the  two-scalar 
formulation  should  therefore  be  a  much  better  approximation  to  the  real 
chemistry  than  in  any  hydrocarbon  flame. 


2.  EXPERIMENTAL 

A  brief  experimental  outline  is  given  here,  with  full  details  elsewhere 
(Marren  et  al.,  1995).  The  methane  fuel  jet  is  diluted  with  air  (77  percent  by 
volume,  this  results  in  a  stoichiometric  mixture  fraction  of  0.29)  and  issues 
vertically  from  a  nozzle  of  diameter  6.1  mm  into  a  filtered,  vertical, 
coflowing  7.0  m/s  air  stream,  at  Reynolds  numbers  from  10,100  to  31,000. 
The  C2H2/H2/air  stoichiometric,  premixed  pilot  flame  flow  rate  is  adjusted 
so  as  to  produce  a  steadily  burning  main  flame  without  significant  local 
extinction.  The  same  fuel  nozzle  is  used  for  the  hydrogen  flame,  but  without 
pilot;  the  Reynolds  number  is  13,500.  To  avoid  overheating  of  the  optics,  the 
H2  flame  is  turned  down  between  successive  laser  shots,  with  a  ‘duty  cycle’ 
of  some  10%. 

Pulses  from  a  flash-lamp  pumped  dye  laser  (3  ps,  532  nm)  are  used  to 
obtain  both  Rayleigh  and  Raman  images.  The  horizontal  beam  is  formed 
into  a  10  mm  high  sheet  with  a  waist  t  =  0.56  mm,  and  aligned  to  include  the 
flame  axis.  The  laser  output  coupler  is  replaced  by  a  cylindrical  mirror 
outside  the  laser  housing  so  that  the  flame  can  be  placed  inside  the  lasing 
cavity,  yielding  an  increase  of  approximately  five  times  in  sheet  energy,  to 
between  4  and  8  J  per  pulse  in  this  experiment. 

The  Raman  signal  is  collected  broadband  with  an  intensified  2000^  pixel 
CCD  camera,  gated  at  <3  ps.  In  the  Rayleigh  images  (collected  with  an 
intensified  576  x  384  pixel  CCD  camera)  the  object  field  is  so  placed  that 
coflow  air  is  included  to  one  side,  enabling  corrections  for  shot-to-shot 
variation  in  beam  intensity  and  spatial  variation  across  the  laser  sheet. 
Images  are  matched,  corrected  for  background  and  throughput,  and 
calibrated  by  the  use  of  exposures  in  fluid  of  known,  uniform  concentration. 
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The  Raman  data  are  smoothed  using  the  contour-aligned  method  of 
Starner  et  al.  (1995). 

When  dealing  with  hydrocarbon  flames,  an  inherent  problem  following 
from  the  single-step  chemistry  assumption  in  the  two-scalar  method  is  the 
loss  of  parent  fuel  to  intermediate  species.  However,  compensation  for  most 
of  the  error  arising  from  this  source  can  be  made  in  the  data  reduction:  the 
difference  between  the  fuel  mass  fraction  in  the  one-step  scheme  and  that 
obtained  from  laminar  flame  calculations  of  strain  rate  representative  of  the 
turbulent  methane  flames  (here  200/s)  can  be  used  to  quantify  the  loss.  In 
this  work  it  is  found  that  the  simple  expression  Rm  -  Rmo(  I  -t-  0.8  i>^)  gives  a 
close  fit  that  is  fairly  insensitive  to  the  strain  rate  over  the  range  of  interest. 
Here,  $  =  (r-rambient)  /  (T’ma.'c-T’ambient),  and  Rmo  and  Rm  are  the  raw  and 
compensated  Raman  intensities,  respectively.  When  tested  in  a  laminar 
flame  of  the  same  composition  (Fig.  1),  it  is  seen  that  most  of  the  ‘shelf’  in 
the  raw  data  is  filled  in  by  the  compensation.  The  level  of  compensation  was 
not  varied  with  the  jet  velocity  as  estimates  indicate  that  the  effect  would  not 
be  significant  for  most  of  the  results  of  interest. 

It  is  common  to  quote  the  sheet  waist  thickness  t  as  the  limit  to  the  scale 
of  measurable  scalar  gradients.  However,  this  poor  resolution  applies  only 
when  the  local  scalar  gradient  is  orthogonal  to  the  plane  of  the  laser  beam 
sheet.  At  the  other  extreme,  i.e.,  when  the  scalar  gradient  is  in  the  plane  of 


FIGURE  I  Compensation  for  fuel  loss 
composition  as  for  the  turbulent  flames. 


to  intermediates  in  a  laminar  methane  flame  of  same 
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the  laser  sheet,  the  resolution  is  limited  by  the  pixel  spacing  s  (here  0.0625 
mm  in  the  processed  images);  this  would  be  the  case  in  a  laminar  jet  flame 
wheii  the  sheet  includes  the  flame  axis.  With  increasing  Reynolds  number, 
the  angle  between  the  scalar  gradient  and  the  laser  beam  sheet  departs  from 
zero,  and  a  worst  case  would  exist  at  very  high  Reynolds  number  when  this 
angle  is  uniformly  distributed  in  all  directions.  For  such  a  case,  by 
integrating  over  all  orientations  a  simple  expression  is  readily  derived  in 
terms  of  the  beam  thickness  t  and  the  pixel  spacing  s  for  the  ensemble 
average  of  the  ‘resolution’  for  the  tangential  direction,  /  for  the 

orthogonal  direction):  (</>)  =  r/2  +  tt  5/4.  Muss  et  al  (1994)  show  that  for 
low  to  moderate  Reynolds  number,  the  conditions  are  much  more 
favourable.  It  should  also  be  recognised  that  no  matter  how  high  the 
Reynolds  number,  in  the  round  jet  flame  the  mean  scalar  gradient  remains  in 
the  plane  of  the  laser  beam  sheet,  i.e.,  the  favourable  bias  will  persist.  On 
this  basis,  we  estimate  that  (0)  <  0.3  mm  in  the  present  work. 

Following  the  arguments  in  Starner  et  al.  (1994a),  for  the  H2  flame  at 
Re  =  13,500  the  Kolmogorov  scale  k  is  estimated  as  0.074  mm  on  the  axis, 
so  {0)//<c  =  4.1.  In  isotropic,  isothermal  turbulence  it  has  been  suggested 
(Tennekes,  1973)  that  energy  dissipation  peaks  at  around  30 This  is 
unlikely  to  apply  to  scalar  dissipation  in  this  non-isotropic,  reacting  jet  flow: 
a  value  of  20  k  for  the  scalar  dissipation  peak  is  indicated  in  Starner  et  al. 
(1994a),  and  resolution  to  5  k  as  sufficient  to  capture  essentially  all  of  the 
scalar  dissipation.  At  the  radius  where  the  mean  temperature  peaks,  k  = 
0.142,  mm  yielding  (</>)//c  =  2,l.  For  the  air  diluted  methane  flame  these 
authors  give  an  estimate  of  0.095  mm  for  the  Kolmogorov  scale  on  the 
centreline  at  jet  Reynolds  number  of  20,600.  The  resolution  in  these  flames 
should  thus  be  adequate  even  at  the  highest  Reynolds  number  of  31,100. 

To  test  the  sensitivity  of  the  computed  scalar  dissipation  to  the  effective 
resolution  of  the  optics  system,  it  is  useful  to  reduce  the  resolution 
artificially  by  increasing  the  pixel  size:  using  the  hydrogen  flame  results,  the 
corrected  Rayleigh  and  Raman  images  are  average  (‘binned’);  then  the  mean 
scalar  dissipation  in  the  radial  direction,  (  Xr),  is  computed.  Figure  2  shows 
that  the  relative  reduction  in  (  Xr)  is  fairly  large  for  the  first  increase  in  pixel 
size,  from  0,0625  to  0.125  mm.  At  this  scale  the  contribution  to  (xr)  should 
be  negligible,  so  this  reduction  is  likely  to  be  due  to  removal  of  pixel-to-pixel 
signal  noise;  it  is  seen  that  the  attenuation  of  (xr)  is  much  the  same  on  the 
axis  and  in  the  hot  zone,  as  may  be  expected  if  it  is  random  noise  only  that  is 
removed.  By  contrast,  on  again  doubling  the  pixel  size  to  0.25  mm,  the 
attenuation  is  much  smaller.  The  relative  attenuation  is  now  also  larger  on 
the  axis,  where  the  integral  scale  is  smaller  (Kelman  et  al.,  1995).  This  is  a 
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FIGURE  2  Effect  of  resolution  (pixel  size  in  the  object  plane)  on  computed  radial  component 
of  mean  scalar  dissipation,  (xr)  in  the  H2  flame.  The  radius  R  is  normalised  by  the  mixture 
fraction  half-radius,  R^.  Average  over  25  images. 


good  indication  of  the  start  of  genuine  attenuation  of  (xr)-  Further  pixel  size 
doubling  to  0.5  mm  results  in  more  substantial  attenuation,  and  again  the 
relative  attenuation  is  largest  on  the  axis  (^5^20  %).  Extrapolation  of  the  data 
at  0.125,  0.25  and  0.5  mm  pixel  size  to  zero  indicates  that  attenuation  at 
0.125  mm  is  less  than  2%  in  the  hot  zone,  and  around  10%  on  the  axis.  The 
results  presented  below  are  given  for  0.0625  mm  size,  and  it  thus  appears 
that  residual  noise  more  than  compensates  for  the  attenuation,  (i.e., 
measured  (x)  is  above  the  true  level)  but  errors  due  to  noise  and  limited 
resolution  should  be  no  larger  than  15%. 

As  an  alternative  measure  of  the  noise  level  on  the  computed  scalar 
dissipation,  the  axial  component  of  x  has  been  measured  in  a  laminar  H2 
flame,  where  the  mean  gradient  of  i  is  negligible.  In  the  turbulent  flame,  this 
component  is  found  to  be  typically  seven  times  the  measured  value  in  the 
laminar  flame.  This  is  consistent  with  the  deduction  made  above  that  the 
initial  attenuation  of  (xr)  in  Figure  2  (pixel  size  increased  from  0.0625  to 
0.125  mm),  around  15  %,  is  due  to  high-frequency  noise  removal. 

Typically,  100  images  are  used  to  calculate  statistical  quantities  in  this 
paper  using  data  from  a  window  1  mm  wide  in  radius  and  8  mm  in  axial 
length  to  give  a  total  of  about  200,000  measurements  for  each  spatial 
“point”.  Errors  in  nongradient  statistics  such  as  means  and  standard 
deviations  of  mixture  fraction,  fuel  mass  fraction  and  temperature  arise 
from  fundamental  uncertainties  in  the  technique  rather  than  on  sampling 
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statistics.  For  means  these  are  estimated  as  90%  confidence  intervals  of 
±0.02  for  fuel  mass  fraction  and  mixture  fraction  and  ±5%  of  the  absolute 
temperature  for  temperature.  For  standard  deviations  the  errors  are  likely  to 
be  somewhat  lower. 


2.1.  The  Two-Scalar  Method 

Full  details  of  the  two-scalar  technique  for  deriving  mixture  fraction  data 
can  be  found  in  Bilger,  (1993c)  and  Starner  et  aL  (1992,  1994a).  Briefly,  in  a 
two-stream,  adiabatic,  reacting  flow,  using  assumptions  of  homogeneous 
inlet  conditions  and  equal  Lewis  numbers,  and  simplifying  the  kinetics  to  a 
one-step  reaction  between  fuel  and  oxidant,  the  fuel  mass  fraction  7/ (here 
defined  as  including  any  dilutant  such  as  70  %  air  in  the  methane  flame)  and 
enthalpy  H=  CpTjQ  are  used  to  form  the  conserved  scalar  (3  =  Yf+  Cp  TjQ, 
where  Cp  is  the  specific  heat  at  constant  pressure,  and  Q  the  lower  heat 
of  combustion  of  the  diluted  fuel.  The  mixture  fraction  ^  is  then  defined  by 

/?-fe  _  Yf  +  Cp[T-T2)lQ  ' 

^  0^-02  Yf,i+c,{Ti-T2)/Q 


where  subscripts  1  and  2  denote  the  fuel  and  air  streams. 

Rayleigh  measurements  are  used  to  obtain  temperature:  cr  /  Ry,  where 

Ry  is  the  Rayleigh  signal.  The  Rayleigh  cross  section,  cr,  the  mole  weight,  W, 
and  Cp  vary  with  ^  and  T,  and  hence  require  an  iterative  data  reduction 
scheme;  predictions  for  an  opposed  flow  laminar  flame  are  used  to 
determine  a  and  W  for  the  burnt  state.  Values  of  Cp  are  based  on  known 
concentrations  of  reactants  and  products.  The  Raman  signal,  Rnif  measures 
the  fuel  concentration.  With  some  arithmetic,  Eq.  (1)  can  be  cast  in  a  form 
which  expresses  the  mixture  fraction  as  a  function  of  Rm  and  Ry: 


(2) 


where  C\  is  a  calibration  constant. 


3,  RESULTS  AND  DISCUSSION 

Computed  profiles  of  mean  and  rms  of  mixture  fraction,  fuel  mass  fraction 
and  temperature  in  the  methane  flames  are  shown  in  Figures  3  and  4  (Favre 
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FIGURE  3  CH4  flames;  radial  profiles  of  mean  and  rms  fluctuation  of  mixture  fraction  ^  and 
fuel  mass  fraction  Yf, - Re=10,100; Re=14,300;- — Re  =  21,300; Re  =  31,100- 


averaged,  denoted  by  tilde  overbars),  at  an  axial  location  of  25  jet  diameters 
(£>)  from  the  nozzle.  The  profiles  of  (|)  and  (f/)are  similar  to  those  for  the 
diluted  CH4  flame  in  Starner  et  al.  (1994a)  but  centreline  values  are  higher, 
and  turbulence  levels  are  lower.  It  is  thought  that  the  main  reason  is  the 
absence  of  local  extinction  in  the  present  case.  (An  extinguished  flame  would 
approximate  to  a  cold  jet  which  would  mix  down  more  rapidly).  The  mean 
temperature  profiles  in  Figure  4  are  quite  smooth,  and  the  peak  temperature 
decreases  as  the  jet  Reynolds  number  increases.  However,  the  rms 
fluctuation  shows  a  local  minimum  at  the  radius  of  mean  stoichiometric 
composition.  This  feature  is  apparent  at  all  Reynolds  numbers  and  is  no 
artefact:  it  can  also  be  seen  in  Drake  et  al.  (1984)  in  the  upstream  region  of  a 
H2  flame  with  jet  Reynolds  number  8,500.  In  downstream  locations  it  is  not 
present,  and  we  may  speculate  that  it  is  the  effect  of  suppression  of  the 
turbulencfe  Reynolds  number  in  the  zone  of  mean  stoichiometric  composi¬ 
tion,  just  inside  the  radius  where  the  onset  of  intermittency  markedly 
increases  the  temperature  fluctuations.  In  the  flame  of  Starner  et  al.  (1994a) 
the  temperature  rms  does  not  show  a  local  minimum;  there  are  two 
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FIGURE  4  CH4  lame:  radial  profiles  of  mean  and  rms  fluctuation  of  temperature,  T. 
- Re=10,100;  ••••  Re- 14,300;  —  Re-21,300; - Re-31,100. 

differences  between  that  flame  and  the  present  ones  which  may  account  for 
this:  firstly,  the  larger  dilution  in  Starner  et  al.  (1994a)  results  in  larger 
stoichiometric  mixture  fraction  which  shortens  the  flame,  so  that  although 
measurements  are  made  at  the  same  axial  position,  it  is  equivalent  to  a 
position  further  downstream  in  the  present  flame.  Secondly,  the  flame  in 
Starner  et  al.  (1994a)  is  half  extinguished  at  the  neck,  which  can  be  expected 
to  alter  the  turbulence  distribution. 

For  the  hydrogen  flame,  corresponding  profiles  of  mean  and  rms  of 
mixture  fraction,  fuel  mass  fraction  and  temperature  are  shown  in  Figure  5, 
at  20  D  from  the  nozzle.  Measurements  in  a  laminar  H2  flame  during  these 
experiments  (not  shown  here)  yield  a  radial  profile  of  ^  that  is  smooth  and 
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^  flame:  radial  profiles  of  density-weighted  mean  ( — )  and  rms  fluctuation  (”•  •  •) 
of  (a)  mixture  fraction  ^  and  fuel  mass  fraction  Yf  \  (b)  temperature,  T. 


has  no  marked  irregularity  such  as  occurs  in  the  methane  of  Figure  1 
around  stoichiometric.  This  is  an  indication  that  the  calculation  of  ^  is  not 
seriously  affected  by  the  limitations  of  the  two-scalar  method.  The  mean 
temperature  profile  in  Figure  5  has  a  marked  peak,  and  its  rms  fluctuation 
rises  steeply  at  the  radius  of  mean  stoichiometric  composition,  as  also  found 
by  Drake  et  aL  (1984).  The  rapid  rise  in  the  rms  of  Tat  RjR^^  1.3  (where 
R^  is  the  mixture  fraction  half-radius)  is  mainly  an  effect  of  intermittency, 
but  may  also  in  part  be  associated  with  a  marked  departure  from 
equilibrium  which  manifests  itself  around  stoichiometric  composition  and 
results  in  a  wider  range  of  T, 

For  both  fuel  types  there  is  evidence  in  the  scatter  plots  of  temperature 
in  Figures  6  and  7  of  marked  departure  from  equilibrium,  as  also  seen  in 
Barlow  and  Carter  (1994).  For  the  hydrogen  flame  the  scatter  plots 
correspond  fairly  well  with  the  prediction  shown  for  a  laminar  counter-flow 
flame  at  strain  rate  1000  s~\  Such  differences  as  can  be  seen  are  thought 
to  be  the  combined  effects  of  higher  (and  fluctuating)  strain  rate  in  the 
turbulent  flame,  calibration  uncertainty  and  the  assumptions  of  one-step 
kinetics  and  equal  species  diffusivities.  The  scatter  plots  for  the  instanta¬ 
neous  scalar  dissipation  that  are  shown  here  assume  that  the  azimuthal 
component  is  equal  to  the  radial  component.  While  this  assumption  should 
be  accurate  for  conditional  and  unconditional  values,  it  will  overestimate 
departures  from  such  means.  This  assumption  is  also  used  in  later  figures 
and  those  results  should  also  be  treated  with  caution. 

The  probability  density  functions  (pdfs)  of  Tin  Figures  8  and  11  show 
largely  the  expected  distribution  in  the  outer  shear  layer,  with  an  inter- 
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‘  Mixture  fraction 

FIGURE  6  CH4  flame:  Scatter  plots  of  temperature,  7,  and  fuel  mass  fraction,  T/,  and  scalar 
dissipation,  x  conditional  on  mixture  fraction,  with  points  collected  equally  from  all  radii. 
Re  =  21,300.  Lines  denote  laminar  flame  calculations  at  strains  10  and  200  s“'. 


mittency  peak  at  ambient  temperature.  The  pdf  of  fuel  mass  fraction  in 
Figure  9  is  somewhat  more  bimodal  in  the  outdr  layer  than  those  of  T.  The 
corresponding  pdfs  of  ^  (Figs.  10  and  11)  are  close  to  Gaussian  on  the  axis,  and 
tending  towards  bimodal  at  the  edge  of  the  flame.  There  is  some  indication  in 
the  methane  flames  (Fig.  10)  of  distortion  of  p{0  around  stoichiometric,  in 
particular  at  low  Reynolds  number,  most  likely  due  to  loss  of  parent  fuel  to 
intermediates  and  the  imprecise  attempt  to  compensate  for  this  effect  (Fig.  1). 

Mean  and  rms  values  of  the  axial  and  radial  components  of  scalar 
dissipation  are  shown  in  Figures  12  and  13.  In  the  CH4  flames  the  axial 
component,  {Xx)i  is  somewhat  larger  than  the  radial  component,  (xr)  OR  the 


FIGURE  7  H2  flame:  scatter  plots  of  temperature,  T,  and  fuel  mass  fraction,  Jy-,  conditional 
on  mixture  fraction,  with  points  collected  equally  from  all  radii.  Lines  denote  laminar  flame 
calculations  at  strain  I000s‘“^ 


FIGURE  8  CH4  flames:  probability  density  of  temperature  F (arbitrary  units)  at  various  radii 
(uncalibrated). 
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FIGURE  9  CH4  flames:  probability  density  of  fuel  mass  fraction  Yf„  (arbitrary  units)  at 
various  radii  (uncalibrated). 


axis,  whilst  in  the  shear  layer  (xr)  grows  to  nearly  twice  (x.v)-  Results  for  the 
air-diluted  CH4  flame  of  Stamer  et  al.  (1994a)  show  a  similar,  anisotropic 
structure.  The  rms  of  the  components  of  x  in  Figure  12  is  around  three 
times  the  corresponding  mean;  this  is  somewhat  more  than  reported  in 
Sdrner  et  al.  (1994a).  It  is  most  noteworthy  that  there  is  very  little  variation 
with  the  jet  Reynolds  number.  This  has  also  been  observed  by  Kelman  and 
Masri  (1996)  for  similar  flames.  For  self-similar  flows  it  would  be  expected 
that  the  scalar  dissipation  would  be  proportional  to  the  jet  velocity  divided 
by  the  jet  diameter  at  a  fixed  number  of  jet  diameters  downstream.  Because 
of  the  coflow,  the  flow  here  is  not  self  similar:  the  spreading  rate  increases 
with  jet  velocity  and  the  mixing  profiles  also  change.  This  probably  accounts 
for  the  lack  of  variation  in  the  scalar  dissipation  values,  but  the  result  should 
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FIGURE  10  CH4  flames:  probability  density  of  mixture  fraction  ^  (arbitrary  units)  at  various 
radii  (uncalibrated). 


FIGURE  11  H2  flame:  probability  density  of  temperature  and  mixture  fraction  (arbitrary 
units)  at  various  radii.  The  radius  R  is  normalised  by  the  mixture  fraction  half-radius,  R^. 


be  treated  with  some  caution.  In  the  H2  flame  of  Figure  13  the  shear  layer 
value  of  (xr)  grows  to  only  1.3  (x.r);  this  is  a  more  isotropic  structure  than  in 
the  CH4  flames.  It  is  not  obvious  why  this  should  be  so,  especially  as  the  H2 
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FIGURE  12  CH4  flames:  radial  profiles  of  mean  and  rms  fluctuation  of  scalar  dissipation: 
components  Xx  and  Xr. - Re  =  1 0, 1 00;  •  •  •  •  Re  =  1 4,300; Re  =  2 1 ,300; Re  =  3 1 , 1 00. 


FIGURE  13  H2  flame:  radial  profiles  of  mean  ( — )  and  rms  fluctuation  ( — )  of  scalar 
dissipation  x-  (a)  axial  component;  (b)  radial  component. 


flame  has  a  rather  low  Reynolds  number.  There  are,  however,  two  marked 
differences  between  these  flames:  in  the  H2  flame  the  density  is  lowest  on  the 
centreline,  and  the  stoichiometric  contour  lies  at  the  outer  edge  of  the  shear 
layer,  whilst  in  the  CH4  flames  the  stoichiometric  contour  (and  minimum 
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density)  lie  well  inside  the  shear  layer  (  Starner  et  al.,  1990).  The  rms/mean 
for  each  component  in  the  H2  flame  is  around  3;  this  is  close  to  the  results 
for  the  CH4  flames  in  this  paper. 

The  pdfs  of  In  (x,0  and  In  (Xr)  on  the  axis  (Figs.  14  and  15)  are  more  peaky 
than  the  lognormal  distribution  and  somewhat  negatively  skewed.  In  the 
shear  layer,  the  effect  of  intermittency  can  be  seen.  A  peaky  distribution  is 
also  evident  in  Rayleigh  scattering  measurements  in  H2  flames  by  Muss 
et  ai  (1994),  but  the  sign  of  the  skewness  for  their  undiluted  H2  flame  is 
positive.  By  contrast,  results  by  Namazian  et  al.  (1988)  in  the  upstream 
region  of  a  cold  CH4  jet  (employing  Raman  measurements)  show  a  close  fit 
to  the  lognormal  distribution.  The  observed  negative  skewness  on  the  axis  in 


FIGURE  14  CH4  flames:  probability  density  of  axial  (left)  and  radial  component  of  scalar 

dissipation  x’ - /?  =  0; - R/R^~  1.07.  Also  shown  is  the  lognormal  distribution  with  mean 

and  rms  as  for  R-Q. 


SCALAR  DISSIPATION  IN  DIFFUSION  FLAMES 


157 


FIGURE  15  H2  flame:  probability  density  of  axial  (left)  and  radial  component  of  scalar 

dissipation  x’ - R-0;  —  R/R^  =  1.07.  Also  shown  is  the  lognormal  distribution  with  mean 

and  rms  as  for  R-0. 


the  present  work  may  be  a  feature  of  the  upstream  distribution  of  scales  and 
scalar  dissipation  in  a  combusting  jet:  on  the  axis  the  scales  are  smaller  than 
in  the  hot  zone  surrounding  the  shear  layer,  and  the  occasional  appearance 
on  the  axis  of  hot  fluid  (with  large  scales  and  low  x)  would  create  a  tail  in  the 
pdf  towards  low  x-  It  should  be  noted,  however,  that  as  limited  spatial 
resolution  will  have  the  largest  averaging  effect  on  samples  with  high  x  (high 
scalar  gradients  being  associated  with  smaller  scales),  the  effect  is  to  clip  the 
pdf  at  high  xi  this  would  contribute  to  an  artificial  peakiness. 

Figures  16  and  17  show  conditional  mean  scalar  dissipation,  {x\v)*y 
conditional  on  mixture  fraction,  and  normalised  by  its  local  unconditional 
mean.  Here  77  is  the  independent  sample  space  variable  corresponding  to  the 
dependent  mixture  fracture  variable,  The  profiles  of  (xl^)*  have  no  clear 
resemblance  to  those  in  cold  flows,  such  as  the  skewed,  monomodal  shapes 
in  jets  by  Kailasnath  et  aL  (1993)  or  the  u-shapes  in  the  grid  flow  in  a  mean 
transverse  temperature  gradient  by  Jayesh  and  Warhaft  (1992).  The  upturn 
in  (xl?;)  at  large  rj  could  be  caused  by  the  rare  appearance  of  fluid  of 
simultaneously  large  ^  and  xi  a  strong  indication  of  this  is  reported  by 
Jayesh  and  Warhaft  (1992). 

Departure  from  equilibrium  (i.e.,  temperature  depression  in  the  scatter 
plot)  is  normally  assumed  to  be  caused  by  locally  high  values  of  x*  If  this  is 
the  case,  one  should  expect  a  clear  negative  correlation  between  T  and  x 
(provided  that  the  measurement  resolution  is  adequate).  This  is  indeed  so:  in 
the  H2  flame  of  Figure  18  the  correlation  coefficient  Ct^v  (conditional  on 
mixture  fraction)  has  a  minimum  of  -0.44  slightly  rich  of  stoichiometric. 
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^^^URE  16  CH4  flames:  Mean  scalar  dissipation  components,  conditional  on  mixture 
fraction  ^  taking  the  value  77,  and  normalised  by  their  local  unconditional  means.  Data  collected 
near  the  axis  and  in  the  shear  layer. =  0-0.25;  — 0.75- 1.25. 

Here  x  is  calculated  by  assuming  that  the  radial  and  azimuthal  components 
are  equal,  i.e.,  x  “  ^Xr  (Starner  et  a!,,  1994a)  and  the  caution  given 
earlier  in  using  this  for  non  mean  quantities  applies.  If  we  assume  that  the 
rms  of  X  is  inflated  by  residual  noise,  as  indicated  in  Figure  2,  then  the 
correlation  would  be  straightened  by  this  effect,  perhaps  compensating  for 
the  error  introduced  by  assuming  that  the  radial  and  azimuthal  components 
are  equal  for  fluctuations.  Outside  the  main  reaction  zone  the  correlation  is 
seen  to  weaken.  This  is  to  be  expected  where  the  flow  is  convection 
dominated  and  the  very  low  reaction  rate  (whether  inhibited  or  not)  is 
insignificant.  The  weak  positive  correlation  for  t]  >  0.2  arises  as  a  calculation 
artefact:  because  finite  mixture  fraction  intervals,  6t],  must  be  used  in 
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FIGURE  17  H2  flames:  Mean  scalar  dissipation  components,  conditional  on  mixture  fraction 
f  taking  the  value  rj,  and  normalised  by  their  local  unconditional  means.  Data  collected  near  the 
axis  and  in  the  shear  layer. - RjR^~0-0.25;  •  •  •  ♦  =  0.75- 1.25. 


FIGURE  18  H2  flame:  (a):  Correlation  coefficient  of  temperature  and  scalar  dissipation, 
conditional  on  mixture  fraction.  The  radial  and  azimuthal  components  of  x  assumed 
equal,  so  that  x  -  Xx'^  2xr*  (b):  Radial  variation  of  correlation  coefficients  of  mixture  fraction 
^  and  scalar  dissipation  components  Xx( — and  xX — )• 


computing  CtJt],  and  because  at  large  77  dTIdrj  and  dxldi]  have  the  same 
sign,  a  spurious  ‘gradient  transport’  correlation  arises  which  can  be 
minimised  only  at  the  cost  of  increased  scatter  as  6t]-^0. 

If  the  assumption  is  made  that  x  and  ^  are  independent  variables,  the  joint 
pdf  of  X  and  ^  can  be  modelled  as  the  product  of  two  separate  pdfs, 
permitting  major  computational  simplification.  The  correlation  coefficient 
(shown  in  Fig.  18  separately  for  the  axial  and  radial  components  of  x) 
appears  to  support  this  assumption  in  the  inner  shear  layer.  For  R/R^>  1.2, 
however,  becomes  substantial.  Similar  results  for  a  CH4  flame  can  be 
found  in  Starner  et  al.  (1994a),  and  for  a  cold  jet  in  Namazian  et  ai  (1988). 
It  should  be  noted  that  a  near  zero  correlation  coefficient  is  not  sufficient  to 
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establish  independence  of  x  and  As  seen  in  Figure  17  the  interdependence 
is  non-linear  and  the  conditional  average  scalar  dissipation  differs 
significantly  from  the  unconditional  average  for  much  of  the  range  of  r?  of 
interest. 

One  important  aspect  of  Conditional,  Moment  Closure  modelling  is  that 
the  number  of  spatial  coordinates  often  can  be  reduced.  It  has  been  shown 
(Bilger,  1993b)  that  in  reacting,  constant-density  flows,  conditioning  on  the 
mixture  fraction  renders  mean  scalar  quantities  independent  of  the  cross¬ 
stream  coordinate.  Recent  development  in  CMC  modelling  (Klimenko 
et  al.,  1995)  indicate  that  the  integral  of  (xh)  across  the  flow. 


f 


{x\v)Prj{vy^r  /  /  p^(v)f^dr 
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FIGURE  19  CH4  flames:  Scalar  dissipation,  integrated  across  the  flow,  and  weighted  by  the 
pdf  of  the  mixture  fraction  (Eq,  3). 
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weighted  by  the  mixture  fraction  pdf,  is  an  important  quantity.  This  integral 
has  been  computed  for  the  CH4  flames,  and  is  shown  in  Figure  19  for  the 
full  range  of  Reynolds  numbers.  Once  again  the  lack  of  increase  in  general 
magnitude  with  jet  velocity  is  unexpected  and  may  be  due  to  lack  of 
similarity  arising  from  the  significant  co-flow.  The  peaks  at  high  mixture 
fraction  probably  arise  from  low  number  statistics.  The  shape  at  low 
mixture  fraction  is  noteworthy  -  an  almost  parabolic  rise  is  indicated.  It 
is  possible  to  obtain  this  weighted  average  of  the  scalar  dissipation  by 
integrating  the  pdf  transport  equation.  Measurements  of  the  axial  velocity 
component  are  needed,  however,  but  are  not  available  at  this  time.  This 
could  give  a  useful  check  on  the  measurements,  as  measurements  of  pdfs  and 
mean  velocity  are  not  sensitive  to  the  degree  of  resolution  of  the  fine  scales, 
as  is  needed  for  direct  measurement  of  the  scalar  dissipation. 


4.  CONCLUDING  REMARKS 

The  results  demonstrate  the  viability  of  the  two-scalar  method  as  applied  to 
hydrogen  and  air-diluted  methane  diffusion  flames.  The  high-energy 
intracavity  design,  combined  with  contour  aligned  signal  smoothing  and  a 
spatial  resolution  of  two  to  four  times  the  Kolmogorov  scale,  permits  2D 
scalar  dissipation  imaging  to  be  performed  with  less  than  15  percent  error 
due  to  noise  and  limited  spatial  resolution.  The  good  agreement  with  point 
measurements  by  others  appears  to  justify  the  assumption  of  one-step 
chemistry. 

The  probability  density  of  scalar  dissipation,  /^(Inx),  differs  from  that  in 
isothermal  flows.  This  appears  to  be  an  effect  associated  with  combustion, 
which  in  this  moderate  Reynolds  number  flame  reduces  the  local  turbulence 
Reynolds  number  and  increases  the  integral  length  scale  around  stoichio¬ 
metric.  More  comprehensive  measurements  will  be  needed  to  resolve  this 
issue. 

The  finding  that  the  scalar  dissipation  does  not  increase  with  jet  velocity  is 
an  unexpected  result  that  must  be  treated  with  caution.  It  is  unlikely  that  it 
arises  from  insufficient  resolution  as  the  measurements  are  too  well  resolved 
to  give  so  significant  an  error.  It  is  more  likely  that  the  result  is  due  to  lack  of 
self  similarity  in  the  flow  and  the  greater  spread  of  the  jet  at  high  velocity. 

Measurements  of  scalar  dissipation  conditional  on  the  mixture  fraction 
show  behaviour  which  is  inconsistent  with  the  assumption  of  statistical 
independence  that  is  sometimes  assumed.  In  the  outer  part  of  the  flow  the 
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conditional  scalar  dissipation  increases  strongly  with  mixture  fraction  at  low 
values  of  mixture  fraction  and  this  is  consistent  with  the  strongly  positive 
correlation  found  between  the  fluctuations  in  dissipation  and  mixture 
fraction.  Near  the  centreline  the  correlation  is  not  strong  but  there  is  a 
significant  nonlinear  interdependence  of  x  and  Results  are  also  presented 
for  the  weighted  integrals  of  the  conditional  dissipation  across  the  flow  as 
are  needed  by  the  conditional  moment  closure  method. 
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